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decay form factors at order p® of chiral perturbation theory which is the next- 
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Leutwyler. At order p® the chiral expansion contains 1- and 2-loop diagrams 
which are discussed in detail. The irreducible 2-loop graphs of the sunset topol¬ 
ogy are calculated numerically. In addition, the chiral Lagrangian produces 
direct couplings with the VT-bosons. Due to these unknown couplings, one can 
always add linear terms in to the predictions of the form factor /_(q^). For 
the form factor /+(g^), this ambiguity involves even quadratic terms. Making 
use of the fact that the pion electromagnetic form factor involves the same q'^ 
counter term, the g'^-ambiguity can be resolved. Apart from the possibility of 
adding an arbitrary linear term in q^ our calculation shows that chiral pertur¬ 
bation theory converges very well in this application, as the O(p®)-corrections 
are small. Comparing the predictions of chiral perturbation theory with the re¬ 
cent CPLEAR data, it is seen that the experimental form factor f+{q^) is well 
described by a linear fit, but that the slope A+ is smaller by about 2 standard 
deviations than the 0{p‘^) prediction. The unavoidable counter term of the 
0 (p^) corrections allows to bring the predictions of chiral perturbation theory 
into perfect agreement with experiment. 


^Now at Fraunhofer Gesellschaft, Institute for Algorithms and Scientific Comput¬ 
ing (SCAI), Schloss Birlinghoven, D-53754 Sankt Augustin, Germany, E-mail: pe¬ 
ter.post@scai.fraunhofer.de 

^Supported in part by Bundesministerium fiir Bildung und Forschung, Bonn, Germany, 
under Contract 05 HT9UMB 4, E-mail: karl.schilcher@uni-mainz.de 


1 




1 Introduction 


The hadronic matrix elements of weak decays constitute a decisive testing 
ground for our understanding of low energy strong interactions. In this re¬ 
spect, the semileptonic if; 3 -decay is one of the cleanest and most interesting 
processes. In particular, it has been stressed that this decay constitutes 
the best source for the extraction of the CKM matrix element |Tus|- On the 
experimental side there exist a number of old high statistics results (some con¬ 
tradictory), and new precision experiments are in progress or already published 
[|| . On the theoretical side chiral perturbation theory Q has established itself 
as a powerful effective theory of low energy strong interactions. Based on the 
symmetry of the underlying QCD, chiral perturbation theory produces a sys¬ 
tematic low-energy expansion of the observables in this regime. Unfortunately, 
because of the non-renormalizability of the effective theory, higher powers in the 
energy expansion require higher loop Feynman integrals and as input an ever in¬ 
creasing number of renormalization constants. The p^-Lagrangian involves ten 
free parameters which were determined in the fundamental papers of Gasser 
and Leutwyler |^. In a more recent analysis the number of independent 
counter terms was reduced to 90. Unfortunately, there is no simple prescription 
how to translate from one base into the other. The series of calculations on 
form factors taken up by us begun before this paper was published. In order 
to preserve the continuity with our previous papers on order p® chiral pertur¬ 
bation theory, we stick to the convention of j^. There is no harm in keeping 
more than the minimal number of operators, as it is, for all practical purposes, 
impossible to determine by the necessary independent experiments such a large 
number of counterterms (including their finite parts), be it 143 or 90. Relations 
between individual physical processes are manifest in any case, as demonstrated 
in iQ and in the discussion below, where a relation of a part of the Ki^ form 
factor to the pion electromagnetic form factor is established. The RTja-decay 
amplitude has been calculated some time ago to [|| and more recently to 

0{p‘^,{md — mu)p^,e^p^) fl. To 0{p^), the contribution of the double chiral 
logs has been calculated in ll0|. 

We present here the results of a full p® and two-loop analysis of the semi¬ 
leptonic Ki 3 form factors, relying heavily on recent progress in the calculation 
of massive two-loop integrals 00- Complementary calculations to this order 
involve the vector two-point function ||l^ and the K 14 form factors |0 . 

As the relevant part of the p®-Lagrangian contains so many unknown pa¬ 
rameters, one may question the usefulness of such calculations. Here are a few 
arguments in favour: 

• In a given class of experiments, such as the electromagnetic and weak 
form factors of the light mesons, only a limited number of renormalization 
constants enter and relations between amplitudes can be tested Q| . 

• The unknown constants enter only polynomially, and precision experi¬ 
ments could separate the unambiguous predictions. 
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• Knowledge of the exact low-energy functional form of an amplitude may 
be important for the experimental extraction of low-energy parameters 
such as charge radii. 

• The results may be used in model calculations which predict the polyno¬ 
mial terms. These calculations can then be compared with experiment. 

• The question of convergence of the chiral perturbation theory may be 
addressed. 

2 The matrix element 


K-mesons can decay into a pion and a lepton pair via the following channels 

{P2) (pi) ( 1 ) 

K°{pi) ■n~ {p2) {pi) viijpv) (2) 

and their charge conjugate modes. The symbol ^ stands for e or We work 
in the isospin symmetry limit (m„ = rrid) where all the hadronic Ki^-decay 
matrix elements are equal. We therefore restrict the discussion in the following 
to Kjg-decay. 

In the Standard Model only the vector current = u'y^s contributes to 
Kia-decay, and the hadronic matrix element has therefore the general form 

(7r“(p2) I I K°{pi)) = {pi +P2)^ f+iq^) + {pi-P2)^, f-iq^) (3) 

where q = pi — P 2 - The dependence of the form factors is usually approxi¬ 
mated by 

/±(a = /±(0)-(l + A±^) . (4) 

The experimental method for the determination of A± consists in comparing 
the measured q^ distribution with a simulation using a constant form factor 
(A± = 0). This approximation could possibly be too crude for future accurate 
data. 

The slope A+ has been remeasured recently in the CPLEAR experiment 
with the result 


A+ = 0.0245 ± 0.0012,t,t ± 0.0022,^,,*. (5) 

This value differs by almost two standard deviations from the previous world 
average [|^ of 

A+= 0.0300 ± 0.0026 (6) 

(which is based on old data of the seventies) and from the prediction of order 
p‘^ chiral perturbation theory Q . 
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The slope A_ can only be measured in if^a-decay, and its status is even 
more controversial. It is common to consider also the so-called scalar form 
factor (because it specifies the S-wave projection of the crossed channel matrix 
element) 

foiq^) = f+{q^) + — ^f-iq^)- (7) 


3 The Lagrangian of chiral perturbation theory 


In the usual formulation of chiral perturbation theory the pseudoscalar fields 
are collected in a unitary 3x3 matrix 

U{x) = exp[i^^] (8) 

where F absorbs the dimensional dependence of the helds, and, in the chiral 
limit, is equal to the pion decay constant, F = 92.4 MeV. The 3x3 Matrix $ 
is given by 


$ = Ao<)ia 


^ 7r°-I-•s/2 7r+ \/2 K'^ \ 

V^7r“ -7r°-I -\/2 

V2A- } 


(9) 


where Xa are the Gell-Mann matrices. 

An explicit breaking of chiral symmetry is introduced via the mass matrix 


/ ml 0 0 \ 

X = 0 0 (10) 

\ 0 0 2m|. - ml J 


where m^r and mx are the unrenormalized masses of the tt and AT-mesons. The 
mass of the ry-meson is given to this order by the Gell-Mann-Okubo relation 


2 ^2 7 2 

"ir, = 3^*- - 


( 11 ) 


The mass term is related to the quark masses by x = const - diag(m„, Wd, TOs) 
with = md- To calculate form factors, we have to include the interaction 
with external boson helds. This is done by introducing gauge helds 


8 8 
= Y^TJ; and ■ (12) 

a=l a=l 

{Ta = Aa/2) with their held tensors 

Ffj,i/ - fl/iG 

— dfj,T 1/ di/T fj, t 
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(13) 

(14) 



and replacing the usual derivative by a covariant one, 

d^U ^ D^U = d^U + iUl^ - ir^U . (15) 

In this way we have extended the global chiral SU{3) x SU{3) to a local sym¬ 
metry. In case of the weak interaction, the external boson is the W and is given 

by 


with 


I, 


•\/2 sin 6 \y 


(W+T+ h.c.) 


0 


/ 0 Vud Vus \ 

r = 0 0 0 

\ 0 0 0 / 


(16) 

(17) 


(18) 


With the building blocks U, and y we can construct the 

Lagrangian of chiral perturbation theory 


C = £(2) + £(4) £(6) ^ ^ ^ 


(19) 


where £(^”1 denotes the most general expression with 2n powers of mass or 
covariant derivatives that is consistent with the symmetries of QCD. For the 
two lowest orders the result is [Q 

£(") = Tr(Z?^I/D'^C/t) + Tr(x[/t + t/yt) (20) 

£( 4 ) = L^{TT{D^UD^^U^)Y+ L2T:t{D^UD^U'^)Ty{D>^UD'^U'^) ( 21 ) 

-h £3 Tr(£»^t/£i'"C/t£)^t7£)‘'I/t) £4 Tr(£»^t7£i^t7^) Tr(xI7^ -h Ux^) 

+L5 Tr(Z?^[/D'^C/t(yC/t + C/yt)) + £g {Tr(x[/t + [/yt)}2 
+£7 {Tr(xt[/ - c/ty )}2 + £g Tr(xt/^xC/t + Ux^Ux^) 

-iLg Tr(£^,Z?^C/£)"C/+ + D‘'U) + £10 Ti{L^,UR'^’^U^) . 


The so-called low energy constants £ 1 ,... ,£10 are unrenormalized coupling 
constants which must be determined by comparison with experiment. The 0{p^) 
Lagrangian was determined in Q. Out of the 143 terms, we reproduce here only 
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( 22 ) 


those which are relevant to semileptonic itT-decays: 

= f3^TT{[D^D,UU{[D^UU[D’'x]^}) 

+ /3i4 Fv{[D^UU ({[^'^X]p, Mm} - {M'^xlm, MpD) 

+ /3i5 {TtHD^UU [D^x]p) Tr(Mm) 

- TtHD^UU Mp) Tr(M'^x]„,)} 

+ /3i6 {TriiD^UU [D^XU Tr([x]p) 

- Tr([i?^C/]„ [X]„) Tr([i^^x]p)} 

+ /3i7 Ty{[D^UU [D^UU Mp Mp) 

+ /3i8 Tt{[D^UU M'^Mm Mp) Tr(Mp) 

+ /3i9 Tt{[D^UU Tr([x]pMp) 

+ P 20 MId^uu Mp) MiD'^uu Mp) 

+ P 21 TtHD^UU Tr([x]p) Tr([x]p) 

+ if322 TtHD^DPUU [[DpUU, M.G'^'lp]) 

+ iP23 Tr{[DpDPUU [[DMU MpG'^'Jp]) 

+ ip24 Tr{[DpUU [DMU {Mp, [G'^'IpD 
+ i/325 MiDpUU Mp [D.UU [G^"]p) 

+ i/326 Tr{[DpUU [DMU [G'^'Jp) Tr([x]p) 

+ 1(327 Tr{[DpUU {[[D.xU [G'^Ip] " [Mm, M.G'^lp])) , 

where we used the notation of § 


Gp. 

— Rpi/U FULpi, 

(23) 

[Mm 

= \{AU^-UA^) 

(24) 

[Mp 

= i(Ac/t + [/M). 

(25) 

with a slight change of notation for the couplings 

Pi = F^Bi (Fearing-Scherer) 

(26) 


so as to use dimensionless quantities. 

The Feynman rules can be derived by expanding U = exp(i$/F’) everywhere 
in £ = + £G) _|_ £(6) identifying the relevant vertex monomials in C. 

Before discussing the Feynman rules in detail in the next section, we would like 
to make a remark on the definition of the form factors. The currents entering 
in i) are defined on the quark level. The connection to the effective theory is 
established by identifying these currents with the Noether currents of the chiral 
symmetry, 

u^pS = , (27) 

where = /(( + r^, a = 1,... , 8, denotes the vector current in the effective 
theory. 
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It is obviously necessary to distinguish also graphically the vertices from 
£( 2 )^ £(4) £(®). We use the conventions: a filled circle • stands for a vertex 

from a filled square ■ for a vertex from and an open square □ for a 
vertex from 

3.1 Pure meson vertices 

If loop corrections to form factor diagrams are considered various vertices enter 
which involve only mesons. In the form factor calculation to 0{p^) five vertices 
from C = will contribute: 

From £( 2 ) we have to consider vertices with 4 and 6 meson fields 

= 2Tr([</),5^</.],^5'^<^)+Tr(x</>") (28) 

and 

1440F4£(2 = 2Ti{d^4>d>^cl)cl)^)-Si:T{d^(t>cl)d>^4>4>^) (29) 

+ 6 Tr(a^((i cj)^) - Tr(x^®) • 

Similarly we have the two- and four-meson vertex from 

= 2L4Tr(a^</)a^</))Tr(x) (30) 

-f 2 L 5 Tr {xd^,(t)d>^(j)) 

- 4L6 Tr(x,^2) Tr(x) 

- 4 L 7 {Tr(x<(>)}' 

- 2X8 {Tr(x'<(<')+Tr(x<(>X</>)}, 

6^"4oI = 6Li{Tr(9,<(.a»}2 (31) 

+ 6 L 2 Tr(a^(/) d,(j)) Tr(a^(/) 5"<()) 

-I- 6 L 3 Tr(a^(/) d^(t) d^fj) 

- 2 L 4 {Tr([0, d,4>\^ Tr(x) + 3Tr(a,0 Tr(x0")} 

- Ls {2Tr(x<(>"5,<^ 5"<(.) + 3Tr(x</> 3^ cf) - Tr(x</> d,4> ^ 

-I- Tr(x d^(j) <t>'^3'"(t)) - Ti'(x (t> d''4> 4>) + 2Tr(x d^4> d''4> 'P^)} 
+ 2Le {Tr(x<(>") Tr(x) + 3[Tr(x<(>')]"} 

+ 8 Lr TrixP^) Tr{xP) 

+ Ls {Tr(x^0'‘) -f 2Tr(x9i>X<(>^) + 3Tr(x«<>^X<(>^) + 2Tr(x^^X<(')}- 

The parameters Li,... , Liq, the so called low energy constants of , are not 
fixed by the symmetries but must be determined by comparing perturbative 
results with experimental data or with models. The low energy constants also 
serve to renormalize the loop diagrams. Therefore they contain divergent pieces 
which, in dimensional regularization, manifest themselves in 4-poles {D = 4 — 

2 e) and which were calculated in Q. 


7 


Finally, the two-meson vertex from enters: 

= -/3i7Tr(a^</)9^0XX) (32) 

-/3i8 Tr(a^(/) 9 ^(()x) Tr(x) 

-/3i 9 Tr(a^^ d^^4>) Tr(xx) 

-/32o Tr(a^^x) Tr(5'"^x) 

-1321 Tr(a^(/) Tr(x) Tr(x) • 


3.2 W-boson-meson vertices 

Every diagram contributing to iFis-decay contains one vertex where the exter¬ 
nal hF-boson couples to the mesons. The Feynman rules of the corresponding 
vertices result from the terms in C that are linear in the gauge fields. Thus, 
the left-handed and right-handed mesonic currents that couple to the external 
pseudo-scalar mesons are given by 

( 33 ) 

r^=F=0 

The result for the currents from reads 


SC 




J. 


R 


=n 


SC 




Tr{T4d^S^,Si]) 


( 34 ) 


48E2 

i 


1440E4 


TrTa{[df,(l>, (j)^] - 5(j)[df,(j>, (j>^](j> + (j>](j>'^) 


where we have only written the terms contributing to iF/ 3 -decay at 0{p^). We 
represent a W-vertex from by a filled circle •. 

The result for the relevant terms of the current from C^'^^ reads 


= 2iLiTT{T4d^S>,S>])Mx) ( 35 ) 

-b iL^ Tr Taix (j) + d^<S x4' - - (S x) 

-iLg d‘'Tr{T4d^4d,S^]) 



for the two-meson-hh vertex and 


12F^ = 24*LiTr(T,[9^<^,<^])Tr(9,0 9"<^) (36) 

+ 24tL2 Tr(T49"0,(^]) Tr(a^(/) O^cj)) 

+ uiLs TiTaiid^cj}, d,(j) 9>}) 

- 2iLi{QTT{Ta[d^(j),(l)])T:r{x(f^) + TrTa{[d^(j),(j)^] - ?>(l)[d^(j),(j)](j)) Tr(x)} 

- iLs {Trra(x[9^(/>,«i'](/>^ + 2x(l>^d^4> (j) - 2 (j)x{df,4>, 

-HTr Ta (2 dfj,(j) x4>^ + 4 (/)X<(' (/) -F 3 (j)]x(ty^) 

-HTr Ta(3 4>'^x[^^l.<i3, <^] + 2 - 4 (() (9^(/) ^x<(') 

-HTr Ta(-2 ((j^x -2.(j)d^(j) (jy^x + ^]x)} 

-h iL95''{TrTa(-3 <i)[d^<j), d^<i)\(j) + d^<i) cj)^ di,(j) - d^4> cj)^ d^(j)) 

+TrTa(2 d^fj)] + 2 [8^4), d^(j)\(j)^ - [d^(j) (j), d„(j) (f>] 

-TrTa{[4df,4,4^0,4])} 

for the four-meson vertex with one tfo-boson. In the diagrams we represent a 
Ifo-vertex from C^^'> by a filled square ■. 

Finally a 2-meson-IF vertex from contributes: 

4F4j^_J/:(6)] ^ 2i(3sTrTa[x,{df.d,<l>,d’'<l>}] (37) 

-f il3u Tr ra(2 X (Z* X - 2 X >9^^ ^ X + XX dyi4' 

- xx4’df,4' +d^4)4)xx-(I)9^4 )xx) 

+ 2i/3i5 {Tr(x0) TTTa[d^4,x] - Tiixd^^) TrT4<^,x]} 

-f z/3i6Tr(x) TTTa[{d^4>, 4>},x] 

+ 2i(3n Tr Ta{4) 9^,4) xx - 9t,4 XX9 + 4>XX 9^i.4> - XX 9t,4 4>) 

-f 2i/3i8 Tr(x) TiTa{4)9f,4) x - 9^,4) x<l^ + 9x9^,(4 - x9t,(44)) 

+ 4i/3i9 {Tr(xx) Tr {Ta[4,9^4]) + 2 Tr(x[x, 4]) Tr(T, 8^4)} 

+ 4i/32oTr(xa^<^)Tr(T4<^,x]) 

+ 4i/32i Tr(x) Tr(x) Tv{Ta[4,8^4]) 

+ 4z/322 TTTa{8^8^[8^8p4,8f’4] - 8p8^[8i,8p4,8f’4]) 

+ 4 z /323 8''8PTTTa{[8p8p4,8^4] - [8^8p4,8p4]) 

- 4i/324 8''Tv{[8p4,8,4]{x, To}) 

-4z/325 8’' Tr Ta{8p4 x9u4 - 9„4 X9p4) 

- 4i/326 8^Tr{Ta[8p4,9A]) Tr(x) 

-f 21827 8''TrTa{[8p4,{8^4,x}] - [9u4A9ti.4,x}]) ■ 

A IF-vertex from L^^'> is represented by an open square □. From the inter¬ 
actions given above, the Feynman rules can be extracted by transforming into 
momentum space and symmetrizing over the meson fields. 
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3.3 Renormalization scheme 


Before evaluating the loop diagrams, we must specify the regularization and 
renormalization scheme. In our calculation we are using dimensional regular¬ 
ization and the so called GL-scheme which is defined in the following way: each 
diagram of order 0{p^^) is multiplied with a factor where = 4 —2e 

is the dimension of space-time and a{e) is given by 

(4^rr(-l + e) = -^, (38) 

£ 

that is 

a(e) = £(l-7-hlog47r)-he^(^-t-i)-)-e>(£^) . (39) 

Because of a(0) = 0 the total 0{p^) result is unchanged in = 4 dimensions. 
The reason for this modification of each diagram is to eliminate the geometric 
factor (47 r)®r(—1 -|- £) appearing in the 1-loop integrals. This renormalization 
scheme is very similar to the well-known MS scheme, where each diagram is 
multiplied by a factor (47r)“^e'>'^ per loop (instead of 6““^'^^). 

The GL-scheme extends the usual 1-loop scheme introduced by Gasser and 
Leutwyler in a natural way. This can be understood by considering the 
renormalization constants Li of : In G-dimensional space-time they have 
dimension D — 4 and their dimension is made manifest by the mass scale p, of 
dimensional regularization: 

L, = . (40) 


Li{p,D) has the same /r-dependence as a 1-loop integral, because Li itself is 
independent of p. It can be expanded in a Laurent series around e = 0 in the 
same way as a 1-loop integral: 


Lt{p,D) = 


L 


(- 1 ) 


+ L^f\p)+eL\^>ip) + Oie^) . 


( 1 ), 


In the usual 1-loop scheme one chooses 

(- 1 ) _ r. 


L 

( 0 ) 


327r2 


L^ip) = Lr(M)-^[l-7 + log47r], 


(41) 


(42) 

(43) 


where T^ are numbers which can be found in Q. The second term in Lf^'^ is 
constructed so that it cancels in the e°-coefficient after multiplication with the 
GT-factor 


Lf^{p,D) = e-“(^)L,(/i,G) 


L 


(- 1 ) 


€ 
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+ L[°^'^\p)+eL\^^’^\p)+Oie^) 

(44) 







with = Lrifi). 


can be treated 


The dimension of the -parameters f3i appearing in 
in the same way: 


Pi = ^ P^i^J.,D) , 


(45) 


where Pi{^,D) behaves like a 2-loop integral. Its Laurent series in the above 
GL-scheme is given by 



3.4 Mass- and wavefunction renormalization 

To order p®, finite S-matrix elements in chiral perturbation theory are obtained 
by multiplying the unrenormalized one-particle irreducible (IPI) Feynman di¬ 


agrams obtained from C = with a factor -s/Z per external 


meson, where Z is the wave function renormalization constant for this meson. 
To be on familiar ground, we start by calculating the mass- and wave function 
renormalization from the renormalized propagator 


t 


(47) 


p2 _ _ y;(p2^ ’ 


where m denotes the bare meson mass and S(p^) the IPI unrenormalized self 
energy which is given perturbatively by 


S(p2) = Si(p2) + E2(p2) + ... . 


(48) 


The leading G(p^)-contribution vanishes, so that Si, S 2 represent the contri¬ 
butions of G(p^) and G(p®). 

From the condition that the renormalized propagator develops a pole with 
residue 1 at where iriph is the physical or pole mass, one derives the 

conditions 


brr? = ml,^ - rrP = S(p2 = mlf,, rrp, F) 


(49) 


Z ^ = 1- E'(p2 = mlf,, rrp, F) 


(50) 


where m stands symbolically for the set of unrenormalized masses and F is the 
unrenormalized pion decay constant (all assumed to be given as functions of the 
renormalized or physical parameters). Perturbatively, Z is therefore given by 


Z = 1 + 6 Zi+6 Z2 + ... 


(51) 


with the and G(p®)-corrections 


SZi = = rrP ,irP, F) 

SZ 2 = Z,2{mpP) + . 
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(52) 

(53) 





where we have expanded around and used the fact that the term 

involving the second derivative of Ei(p^) vanishes, i.e. i® independent of 
. In eq. (^) the unrenormalized quantities rn?,F must be expressed by their 
physical values according to eq. (|^) and eq. below. 

Each external meson propagator must be multiplied with a factor 

^/z = i + ^ + (54) 

We have to calculate for tt-, itl-mesons in order to determine Zt^,Zk 

and (5m 5m]^. The results are well known [||: 

+ 2Ao(m2) (55) 

+ 24L4(2m|^ + ml) + 24L5m^} 

8 Z^ = ~^^{^o(m^) + 2Ao(m^) + 44o(m^) (56) 

+ 32L4(2m|- + m^) + 32L5m|-} 

^ - 3m2 Ao(ml) (57) 

— 48m^(2m|- + m^) L 4 — 

+ 96L6m^(2m^ + m^) + QGLgm^} 

SrriK = Y^{-4wK^o(m2)-96L4m^(2m| + m2) (58) 

— 96L5m^ + 192LQmj^{2mj^ + m^) + 192L8m^} 

The function Ao{m‘^) is the standard tadpole integral 

where Z? = 4 — 2e is the dimension of space-time. 

In addition, we need the renormalization of the pion decay constant 

5F = F^-F , (60) 

where E,r is the physical pion decay constant. We only quote the result Q: 

5F = ■^{Ao{m\) + 2Ao{ml) + 8L4(2m^ -b ml) + SL^ml} . (61) 

It should be noted, that the renormalization constants 5m? and 5F defined 
above are finite. The divergences and scale dependence of the loop integrals are 
canceled by similar factors in the counter terms Li from 

The self energy diagrams contributing up to two loops and to order p® are 
given in hg. 1. External legs are fixed by the mesons considered, while one has 
to sum over all possible internal meson lines. 
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We call a two-loop diagram ’’reducible”, if the two loop integrations decouple, 
i.e. if they are given by a product of one-loop integrals. Otherwise they are 
called ’’irreducible”. The C>(p®)-correction of the wave function renormalization 
Z consists of three parts 

(5^2 = + bZY"'^ + bZ^ (62) 


which are given below for tt and K. 

The one-loop and the reducible two-loop diagrams of fig. 1 yield for the pion 
wave function renormalization Z'^ an O(p®)-contribution 


2880 [2Li+4L2 + T3]^2(m2) (63) 

-h 2880 [4L2 -f La] ^ 2 (m|) -f 960 [ 3 L 2 -f L 3 ] A 2 (m 2 ) 

-h 960m^ [(2m^ -h - 2Lq) + ml{L^ - 2LY\ A 0 . 2 YI) 

-\- 480m^ T 2m^)(L4 — 2Lq) -\- m^(L5 — 2L8)] AQ.2ijn’Y} 

-\- 720771^ [ 12 L 1 -|- 4 L 2 -f 6 L 3 — 6 L 4 — 3 L 5 ] j4o(^t) 

-I- 480 [24m|.Li -|- G-mj^Ls — 12mj^L4 — (m^ -|- 2 to|^)L5] Ao{mY 
+ 80 [2(4m|; - ml){ 6 Li + L 3 - 3L4) - SmlL^] Ao{m'^) 

+ 20 ml [3^o(to^) - -4o(to^)] Ao. 2 irnl) + 20m|- Aoim^) Ao- 2 imY 
- 15 Ao{mY Ao{ml) + 6 y4o(m|-)^ -h 9 Ao{ml) Aoim^} 

+ {szTf , 


where SZf is the 0 {p'^) result from eq. ( p5[) and the additional functions Ao ;2 
and A 2 are related to the tadpole integral Aq and the dimension D of space-time 
(cf. eqs. (1^) and (1^) in appendix 



1 ^ Ao(m^) and A 2 {iY) = 4lo(m^) . 


(64) 


For the irreducible two-loop contributions of fig. 1, which involve higher 
transcendental functions, we only quote the exact result for the divergent part 
and a numerical result for the finite part. The latter involves an arbitrary scale 
/r which cancels in the final answer for the form factor. For the choice 


pL = mp = 770 MeV, mK = 495 MeV, = 0.28 = 92.4 MeV, (65) 

ml being an abbreviation for the Gell-Mann-Okubo term ^m\- — ^ml, and the 
definition 


Lg(m2) = log(^)+7 
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we obtain 




2.015180506- 


78m]^ + 22 m%m^ + 89m^ 


72(47r)4i^4 £.2 

288(47T)4j^4g {^452TOj^ - 36 to^to 2 + 1163m^ 
—4(24m^ + 14 to|^to^ — 5m^) Lg(m^) 
-528m^ Lg(m%) 

-UmlilOm^ + 61ml) Lg(m^)| . 


( 66 ) 


Similarly we obtain for the reducible part of the kaon wave function renor¬ 
malization 

= 3^{4320 [44.2 + 1^3] ^2(^2) (67) 

-k 2880 [4Li -k IOL2 -k 3L3] 2l2(m^) -k 480 [I2L2 + ^3] ^2(^2) 

-|- 720to^ [(2to^ -|- ml){L4 — 2 Lq) + m^(L5 — 27^8)] 7lo;2(m.^) 

-|- 1440m^ [(2771^ -|- m^){L4 — 2 Lq) — 27^8)] Ao-,2{mj^) 

+ 80 [(4to^ - mD^L^ - 48(m^ - mD'^Lj 

-\- 3(4m^ — ml)(2ml( + ml){L4 — 2 Lq) 

- 6{8m% - 8mj^ml + 3ml)Ls] Ao; 2 (to^) 

+ 720 [6m2 (7.3 - 2L4 + 4Li) - (2777^ + 777^)7.5] ^0(^77^) 

-I- 1440777^ [167^1 -I- 4 L 2 + 6 L 3 — 87^4 — 37 ^ 5 ] Tip( 777 ^) 

-k 80 [2(4777^ - ml){5L3 - 6L4 + 12Li) -f 3(2777^ - 7777|:)L5] ^0(777^) 
-k 45777^ ^ 0 ( 777 ^) ^0;2(777^) “ 45777^ ^o(777^) ^0;2(777^) 

-I- 120777|^ Ao{mli) Ao; 2(777^) - 15777^ ^o(777^) ^0;2(777^) 

-k 60777^ ^0(777^) ^0;2(777^) + 5(7777^ - 16777|;) Ao(777^) 2lo;2(777^) 

-k 45 2lo(777^)^ - 18 Ao( 777^)^ - 27 zlo(777^)^ 

- 54 2lo(777^) Ao{ml) - 27 Ao{ml) ^0(777^) -k 81 ^0(777^) ^0(777^)} 

+ ((5^3^) 2 . 


The irreducible diagram yields 


^‘irred 


2.8602989531-k 
1 

288(47r)4i^4g. 


1427771;- -I- bml^ml -I- 42777^ 
72(47r)4F4 £2 

11769777^ -f 186777|-777^ -|- 624777^ 


( 68 ) 


- 4(56777|f — 26ml^ml + 8ml) Lg(777^) 

- 24777 I; ( 38777 I; -I- 3777^) Lg(777|;) 

- 36777^(2777^ -f 9777^) Lg(777^)| . 
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Finally, the contribution arising from the -constants of Eq. are 


^{ 4to^ /3i 7 -h 4m^ (2 to^ -h ml ) Pis (69) 

-|- 4(4to^ — Aml^rnl 3ml) Pw 
+ 4(4to^ -h ^mj^ml + ml) P 21 } 

■^{A{2m\ - 2m\ml + ml) Pn (70) 

+ 4mj^(2mj^ + ml) Pis 
+ A(Am% - Am\ml + 3ml) /3i9 
-I- 4(4 to^ -I- Amj^ml + ml) P 21 } ■ 

If the unrenormalized contributions of the order p^, p'^ and p®-diagrams of fig. 2 
are denoted as Ag/, Ai/, and A 2 /, then, with the mass- and wave function 
renormalizations given above, the renormalized form factors read 

/ = (Aof + Aif + A2f) (71) 

= Ao/+{Ai/ + Ao/(i<5Zf+ i<5Zr)} 

Z Z m2=m2^-<5m2, F = F7r-5F 

+ A 2 f+ Aif{^SZl^ + ^6Z^} 

+ AoJi^SZ^ + + ^SZl^SZ- - i(5Zf )2 - ^(SZ-)^} 

+ o{p^) 

where / stands for f±. 

3.5 £^®)-contributions to the form factors 




5Z^[L^^^] = 


In every order of chiral perturbation theory there appear new operators with 
a priori unknown coefficients. The ten constants of are by now all fixed 
by experiment, but little is known about the 143 constants of Out of the 
latter, only 11 enter in semileptonic F-decay. There are two sources which lead 
to £*^®^-contributions to the form factors: one is the 0 (p^) tree graph of fig. 2, 
and the other one is the 0 (p^) tree graph of fig. 2 with the 0 (p^) wave function 
renormalization. Since to 0{p^) 

Ao/+ = 1 (72) 

Ao/- = 0, 
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(b) 



Figure 1: Self energy diagrams up to order p^. £^^^-vertices are denoted by filled circles (•), 
£(4)-vertices by filled squares (■), and an £(®Fvertex by an open square (□). 
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the total contribution involving -constants (ordered in powers of q^) to the 
form factors /+ and /_ becomes 

^ {(322 + fe} - ^{/322(m| -f ml) - 2(32imj^ 

-(325ml - /326 (2 tok + ml) + (327{m]i + ml)} 

+ ^l 3 ii{m\ - ml)"^ 

~ ml){Ps + P22 + P23} 

2 2 

+ pi "" [Wsimli + ml) - 2(3iQ{2m\ + ml) 

-2(3is{2m\ + ml) + 2(322{mj^ + ml) - 4:(32im\ 

-2(325ml - 2(32&(2m\ + ml) + 2(327{m]i + ml)} 

The q'^-term of A/+[£(®)], i.e. 2q^{(322 + (32 ^)/, is the same as that for the 
electromagnetic form factor of the charged pion (and kaon). One can therefore 
use data on the second derivative of the pion electromagnetic form factor at the 
origin, together with the 0 {p^) loop calculation, to determine the combination 
(322 + (323- Details are given in appendix]^. We find for the e°-part 

{P22 + / 323 )'‘’^’“ = 0.61 X 10-4 (75) 

where we have used the GL-scheme (^) and chosen the mass scale p = 770 MeV. 
To these counter term contributions, we have to add those of the loop diagrams 
involving and £^4)-vertices. 

3.6 Divergent two-loop contributions to f± 

We start with an analysis of the divergent two-loop contributions of fig. 2 to the 
form factors. The 0{p^) pole terms in e have to cancel in the sum of all loops 
and the tree graphs with an £(®^-vertex: 

+ A/±[/:(6)]) = 0 . (76) 

Since the £^®btree graph part A/± is a polynomial in masses and momenta, 
it follows that the loop part must also be polynomial in masses and 

momenta, i.e. it cannot contain any logarithms thereof. This condition offers a 
good check of the calculation. In fact, we find that in the sum of all 0{p^) loop 


(73) 


(74) 

- 4/3i7to|- 
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(3f) 


(4) 


(5a) 


(5b) 


(5c) 






Figure 2: The diagrams for the K 13 form factor up to order p®. £^^Fvertices are denoted by 
filled circles (•), yC^^Fvertices by filled squares (■), and an £(®Fvertex by an open square (□). 
Diagrams (0) to (3f) are referred to as ’’reducible”, diagrams (5a) to (5c) as ’’irreducible”. 
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diagrams any non-polynomial terms in masses and momenta cancel: 


Ap6°^/+ = 95(4^^j^)4g2 -H 21m^ - 28m^g^ - 3g^} (77) 

+ 576(4^)2 (-Wk + - ml + 

-I- 384L^°^(-m|. -I- 2mj^rnl - - ml - Sm'lq'^ -|- q'^) 

+ 64 + 2m]^ml + — ml + Qm^q^ — 9q'^) 

+ 768 Lf\3m% - 6mj^ml + m\q^ 4- 3ml + m^q^) 

+ 576 L5°^(—3to^ -I- 6m\ml + mj^q'^ — 3ml + m^q^) 

+ 288 Lf\^{lm\ + 5m^ -|- q^) 

- [70m^ - IdOml^m^ -I- 167m^g^ -I- 70m^ -|- 125m^g^ - 12g^]} 

(47r)^ 

and 

Ap6°^/- = 48(4jp)4g2 i^Q”^^ + ^^m\ml - 27m\q^ - 47ml + ‘27mlq^} (78) 

+ 570 ( 47 ]-) 277 ' 4 £' l '^®®+ ^mj^q^ + 7ml — Sm^q^) 

+ 384L^°^(-19to;|- -I- 7m\q^ + I9ml - 7m^g^) 

-I- 64L^°^(—145m^ -I- 48m']^ml + 41mj^q^ + 97ml ~ ^^m^q^) 

+ 384 Lf\-26m% + 9m^m^ -I- 17m^) 

-h 64L®(-m^ -t 2m‘^iiml + 27mW -ml- 27mlq^) 

4- 6912 L® (2to^ — mj^ml — ml) 

+ 9216 L®(m^ — 2m\ml + ml) 

+ 2304L®(5to^ - 4TO^m^ - ml) 

+ 288 L^g\—7m% + 2to^to^ — -I- 5ml + m^q^) 

+ - 204m|;TO^ - Slml^g^ -I- 193 to^ -|- 51m^g^]} . 

This is not the case for the group of reducible resp. the group of irreducible dia¬ 
grams alone, only in their sum. are the e°-coefficients of the constants, 
cf. eq. d^. 

The divergent parts A/±[£(®^], given explicitly in terms of the £(®)-constants 
j3j in eqs. (|^) and (^), have to be the negative of these expressions so that 
the whole 0{p^) prediction is finite. Since the £*^®^-constants /3j do not know 
anything about the masses mj, m\, they must appear in ( fz^ ) and ( [T^ ) in such 
a way that the mass dependence of the divergent parts A/±[£(®)] is produced 
from the explicit masses in (|zl)> dzl) alone. In other words: the £*'®^-constants 
themselves cannot contribute any masses. In appendix we list the resulting 
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divergent parts for the relevant £(®)-constants fij. The fact that they are inde¬ 
pendent of the masses, is a consistency check between our calculation and the 
Lagrangian from Q. 


3.7 Reducible loop diagrams 

The reducible diagrams of fig. 2 can be expressed in terms of one-loop integrals. 
In case of two loops they are of the form 


2D f d^l V{k‘^,kpi,kp2,kl,Pjpi,lp2) 

J (27r)^ P 


and in case of one loop 


f d^k V{k‘^,kpi,kp2) 

J P 


(79) 


(80) 


Here H is a polynomial of its arguments and P represents a product of propaga¬ 
tor factors that depends on the topology of the diagram. A first simplification 
of these integrals is achieved by replacing certain factors in the numerator ac¬ 
cording to 


k^ = P{k,m^)+m^ (81) 

kq = ^[P{k + q,m^) + — k^ — q'^] 

f = P{l,m^) + m^ 

= ^[P{1 + +m'^ - - q'^] , 

where P(k,rn?) = k'^ — m? + zO. The remaining integrals can be expressed 
through the 1-loop one-point function Ao(to^) of eq. (|^), the 1-loop two-point 
function 


Bo{q'^,ml,ml) 


A* 


4-D 


d^k _ 1 _ 

z(27r)'° [(Ac -I- qY — mf][fc2 — m|] 


(82) 


and tensors and mass derivatives thereof. 

In a two-loop calculation these functions have to be considered up to order 
(where H = 4 — 2e). The results are given in appendix The reducible 
contributions to /+ are given explicitly for each diagram in appendix 


4 Irreducible two-loop diagrams 

In the irreducible diagrams 5a, 5b, 5c of fig. 2 the two loop integrations are not 
independent of each other as they were in the reducible graphs. That is why 
genuine 2-loop functions enter the stage which cannot be expressed by 1-loop 
integrals only. 
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Inserting the Feynman rules yields integrands with a similar structure as in 


V{P, kpi,kp 2 , kl, f, Ipi, lp 2 ) 


P{k +Pi,to§)P(A: + p 2 ,m\)P{k + l,m\)P{l,m1) ’ 

where is a polynomial of degree equal to the number of vertices, 
celing factors via 


(83) 
After can- 


kpi 

= \[P{k + Pi,rnl)+ml-k'^- pI] 

(84) 

kp 2 

= ]^[P{k + P 2 :rnl)+ml-k'^- pI] 


kl 

= ^[P{kPl,ml) + ml-k^-t^] 



= Pi^rriTl + ml 



we are left with reducible integrals which can be calculated analytically, and with 
some genuine 2-loop integrals of the swnset-topology, i.e. the 3-point functions 

r7~i / 2 2 2 2 2 2 2\ _ 

7ai.a2./3(9 ;Pl,P2;"ro,Wl,?7l2>'^3) = 


8-2D f d^k d^l _ {lp,r- {lp2r- {Pf _ 

^ J z(27r)'° i{27r)^ P{k + pi,ml) P{k + p 2 ,ml) P{k + l,ml) P{l,ml) 

and the 2 -point functions 




8-2D 


d^k d^l 


{lp)°‘ (k^)^ 


J z( 27 r)^ z( 27 r)^ P{k + p,rn\) P(k + l,m\) P{l,m\)' 

( 86 ) 

In diagram (5b) of fig. 2, nine different mass flows of intermediate mesons 
must be regarded: 


A(5'')/± = 


.(5b) 

^K-ktttz 

,(5b) 


/±+ A 


(5b) 


KttKK 


f± + ^k 1 J± 


+ K7.Kf± + A 


-I ^(5b) r 

^ ‘^r]KKriJ± 


(5b) 


n/± + A, 


KTzrjr] J 

(5b) 

tiK-kK 


f± 


A (5'’) f, 

^Kr)KKJ^ 


^C5b) r 


(87) 


( 88 ) 


where means that mesons r and s couple to the IF-boson and the other 

two lines are mesons of type t and u. Each mass flow is handled separately, and 
its contribution to the Kii, form factor is expressed in terms of the basic 1 - and 
2-loop functions A, P, Pa./s, 2 ai,a 2 ./ 3 ) where for the latter at most the tensor 
indices 


<S' 2 . 0 , S'!,!, Pi^o, 5'o.O, 5'o,l, 5*0,2, To. 0 , 0 , To. 0 , 1 , 70.0,2, To,0.3, Pi, 0 . 0 , 
Pi.0,1, Pi,0,2, Pi,1.0, Pi,1,1 


are needed. Except for special kinematic situations the genuine 2-loop inte¬ 
grals 5a ,/3 and Pai,a 2,/3 cannot be calculated analytically. In El we describe 
the method how we calculated them by splitting them up into one part which 
contains the divergence and can be evaluated analytically, and a second part 
which is finite and can be done numerically. 


(85) 


21 









5 Finite contributions of the loops 


After presenting the results of the pole-terms of the loop-diagrams we now 
come to their finite parts which contain the actual physical information. We 
will present the results, which can only be given in numerical form, graphically 
and as interpolation polynomials. We use the GL-scheme discussed above at a 
scale fi = rup = 770 MeV, the masses given in (^, and the following values for 
the finite parts of the -constants: 


i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

104 ^en 

7±5 

12 ±4 

-35± 13 

-3 ±5 

14 ±5 

-2 ±3 

-4 ±2 

9±3 

69 ±7 


(89) 


The -coefficients of Li are not new degrees of freedom, but always appear in 
combination with certain G^-^^-parameters. Therefore, we define 

= TOp) = 0, (90) 

cf. (^^. As our kinematical range we choose momentum transfers < 

< (jriK — iTT'Tr)'^, which is certainly within the range of applicability of chiral 
perturbation theory. The results of the 0{p'^)- and the G(p®)-loop contributions 
(reducible and irreducible) for the form factors /+(g^) and are plotted as 

a function of a; = /m^j^ in figs. 3 and 4. We observe that for /-|_ the reducible 

and the irreducible O(p®)-contributions cancel almost completely. For /_ the 
irreducible loop corrections are very small. For both /+ and /_ the 0{p^) loop 
corrections are essentially linear functions of i.e. they create only small 
nonlinear contributions to the form factors. 

The interpolation polynomials for the loop corrections at order p® of the 
form factors G (i.e. the 0 {p‘^,p^) parts of the full renormalized form factors 
given in eq. (|7l|), but without the tree graphs from read 

A‘°°Pf+ = a:(-0.0101 -k 0.0009 a: -k 0.0054 a:^ -k 0.0007 a:^) -k Ap°°^/+(0) (91) 
= a:(0.0185 -k 0.0001 x -k 0.0022 + 0.0008 x^) + (0) , (92) 

X = q^fmj^ e [0,0.5]. For completeness, we quote the results for /±(0) which 
come from the 0 {p*) and 0 {p^) loops: 


aXV/+(o) 


A loop 


/-(O) 


Ap4/+(0) + A;f■ '°°"/+(0) + A;r‘'- '““^/+(0) (93) 

-0.0229 -k 0.00937 -k 0.00853 = -0.0050 

Ap4/_(0) + '““P/_(0) + '™^/_(0) (94) 


-0.1836 -k 0.0832 - 0.0533 = -0.1537 . 


Here, denotes the contributions of the irreducible diagrams (5a), 

(5b), (5c) from fig. 2 and the irreducible part of the wave function renormaliza¬ 
tion (diagram (h) from fig. 1). A^e"^' stands for the remaining loop contribu¬ 
tions which all come from reducible diagrams. To these results, the contributions 
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Figure 3: Loop contributions to the form factor /+(g^). The solid line is the leading order 
0 {p^) result, the dotted curves denote the corrections at order p® due to the reducible resp. 
irreducible loop diagrams in the GL-scheme (on the one hand diagrams (2a)—(3f) and the 
contributions of the renormalized leading order diagrams (0), (la)-(lc), on the other hand 
diagrams (5a)—(5c), cf. fig. 2). The contribution of diagram (4) is missing: it would add an 
arbitrary parabola through the origin. 


-/-(o) 



Figure 4: Loop contributions to the form factor /_ (g^), in analogy to fig. 3. The contribution 
from the G^^^-tree graph (4) is missing. It would add an arbitrary straight line through the 
origin. 
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of the constants, which occur at 0{p^), must be added. These can either 
be determined by experiment or by model calculations. 

The 0{p^) loop contributions still depend on the mass scale fi. This depen¬ 
dence follows from their divergent parts which are given in eq. ©• For /+(0) 
we find a value of 0.0617 at the scale p = 500 MeV and a value of —0.0364 
at the scale p = 1000 MeV to be compared with the above value —0.0050 at 
p = 770 MeV. It is seen that the loop contribution to /±(0) is rather sensitive 
to the choice of p due to cancellations of the reducible and the irreducible con¬ 
tributions. This /r-dependence is of course canceled by the /r-dependence of the 
g® counter terms. Eq. (|^) and (|^ only serve as an indication of the size of 
the effects to be expected. 

At the present stage the predictive power of our 0{p^) calculation lies only 
in quantifying a deviation from linear form factor rise. From eqs. (^3|) and (H) 
we find a nonlinear contribution to /+ of 

= 0.10 ^ , (95) 

whereas /_ contains negligible nonlinearities in the relevant kinematic range. 
The effect of the quadratic term in /+ is essentially a lowering of the parameter 
A+ defined in (|^). In our fit we find 

A+ = 0.022 (96) 

(cf. fig. 5) as compared to the linear fit A+ = 0.0245 in 0. 

We conclude this section with a short discussion of the errors of the 0{p^) 
correction. Apart from the general ambiguity due to the counter terms, 
errors arise from the constants which appear in the 0 {p^) corrections, and 
from the one constant, associated with the q‘^ term of /+(q^), which can 
be extracted from the pion electromagnetic form factor. Although these errors 
could be as large as 10%, they are irrelevant to our result as the total 0 {p^) 
effect is small, see fig. 5. 

6 Analysis of results and conclusions 

We have calculated the 0{p^) contribution to semileptonic Ki^ decay in SU{3) x 
SU (3) chiral perturbation theory. This is an effective field theory, so that there 
appear new operators with unknown couplings in each order of perturbation the¬ 
ory. For the Ka form factors f±{q^) this means that the constant and the linear 
term in q^ are not determined by the theory. The g'^-counter term, however, 
is the same for the semileptonic form factor f+{q^) and the electromagnetic 
form factor of the charged pion. It can therefore be extracted by comparing 
the 0 {p^) chiral perturbation theory with low energy data on the pion electro¬ 
magnetic form factor. The details are described in appendix ^ There is no 
(/■^-counter term for the form factor f_[q^). 
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Figure 5: The kaon semileptonic form factor /+ to O(p^) and 0(p®) chiral perturbation 
theory versus experiment. The data are from the CPLEAR experiment 0, the upper dotted 
curve is the leading order 0(p^)-result of chiral perturbation theory, the solid curve is the 
full result at order p®, and the lower dotted curve is the linear fit from 0. The slope at the 
origin of the 0(p®) result has been fitted to the data, its second derivative comes from the 
electromagnetic pion form factor. 


We have found that the 0{p^) loop corrections are essentially linear in q^. 
Thus, the only nonlinearity is the g^-contribution to /+ coming from the 
tree graph which is related to the electromagnetic pion form factor. 

It is interesting to consider the reducible and the irreducible 0{p^) loop 
results separately as plotted in figs. 3 and 4. It is clearly seen that for /+ both 
terms cancel almost exactly in the physical range of iF /3 decay. It should be 
kept in mind, however, that an arbitrary linear term can always be added to 
the 0{p^) predictions. For /_ the irreducible contributions are very small. 

The results obtained are interesting from the following points of view: 

• The convergence of chiral perturbation theory in semileptonic i^js-decay 
is established. It turns out that in decay the 0(p^) corrections are 
small. This need not always be the case |^ . 

• For the form factors, the deviation from a linear rise in is small, but 
not negligible (for /+). The result of the nonlinear contribution to /+ is 
effectively a lowering of the parameter A+. 

• Our method of calculating the irreducible two-loop diagrams of high ten- 
sorial rank and involving three different masses can also be applied else¬ 
where. 

• The wave function renormalization constants calculated here can be em¬ 
ployed in other processes. 
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• The divergent parts of the relevant -parameters f3j appear in other 
processes as well and can be used there as a check. 


A 1-Loop integrals 


In this appendix we reproduce the well known one-loop results. With the defi¬ 
nition 


Ttl^ 

Lg(m") = log (^)+7 


( 97 ) 


one has 


^ J i 


d^k 


(27r)'° — w? 


( 98 ) 


, 1 


(47r)2 ^ £ 


{1 + 1 _ Lg(m^) + e[l + ^ - Lg(m2) + + 0{e^)} , 


Lg(m^)^ 


and 


D ^ 2 „ 2 „ 2 ^ 4 -D f 1 

Bo{q ,m^,m2) - M J ^^27r)D [{k + q)^ - mlW - 


( 99 ) 




+ b^°\q‘^;ml,ml) + eb^^\q‘^;ml,ml) + C’(e^), 


with 


2 2 \ 

;mi,m2) = 


b^^\q'^;rnl,ml) = 


1 


(47r)^ 

1 

(47r)^ 

1 


( 100 ) 


|2-Lg(TO?)-f ^Xjlog(l- —)j 

L .7 = 1 ^ 

Lg(mf)2 


.2\ 


( 47 r) 


2(4+ - 2 Lg(mi) -f 


1 T 1 

+ log(l- )[4 - 2 Lg(g 2 _ zO) - log(l - Xj) - \og{-Xj)] 

j=i 

-Xi[log(l - Xi) log(l - X2) - log(-Xi) log(-X2)] 

+{xi - a;2)[log(x2 - Xi) log(l - —) - Li2 (^^——) + Li2(-^^)]} 

X2 Xl — X2 Xl — X2 
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where D = 4 — 2e is the dimension of space-time, X\j 2 are given by 


Xl,2 = 



( 101 ) 


A being the Kallen-function A(a, b,c) = (a — b — c)^ — 46c, and where we have 
expanded to order All masses carry an infinitesimal negative imaginary part. 

Tensor integrals can be reduced to scalar ones by decomposing them with 
respect to Lorentz covariants. The notation is 


Ar(m^) = coefficient of the tensor integral with r momenta fc ( 102 ) 
in the numerator (r even) 

= coefficient of the tensor integral with r momenta/c (103) 
in the numerator, and s factors of on the rhs. 


e.g. 


r d^k kf^k'' 

J i{ 2 'K)^ — w? 




(104) 


and 


A* 


4 -D 


9- 


4 -D 


9 


4 -D 


d^k k^^ 

i{2Tr)^ [{k + q)'^ — ] [k9 — m|] 

d^k fc'" fc'" 

i{27r)^ [(fc -I- 5)2 — ml] [fc^ — m|] 
d^k fc'" fc"" kP 

i{2Tr)^ [(fc -I- ( 7)2 — mf] [fc^ — mf] 


q^Bio 


(105) 


q^q'^B^o + g^’'B2i 

qPq-'qPB^o + {g^’'q^ + g^^q'' + g''>>q>^)B^i. 


A 2 is given in terms of Aq in eq. (| 6 ^, and the functions Bj-s can all be related 
to Bq (and Aq). For the form factor /+ we need 


B2i{q'^;ml,ml) 


B3i{q^\m\,ml) 


4g2(l - -ml- q^)AQ{ml) (106) 

-fi {ml - ml- q^)Ao{ml) 

+ ,rnl,ml)Bo{q^-,ml,ml)} 

8q4{l- -^ 1-^2 + ‘^mjml + imjq^ - ]Ao(to^) 

+ [A(g^mf,TOi) -f ]Ao{ml) 

+ {ml - ml- q^)\{q^, ml, ml)Bo{q'^; ml, ml)ll07) 


In addition, there are 1-loop integrals which involve higher powers of prop¬ 
agators. These are obtained from the formulae above by differentiation with 
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respect to 


Ao-nirn^) = 

n.T 5^1; ^^2) ~ 


1 


in—1 


(n — 1)! ''dm?' 

1 


Ao(m^) 


(108) 


rs-u,u,yH - (^^_1)!(^2_1)! L ^^2 


lit'-'“t”») (1»9) 


B Reducible contributions to the form fac¬ 
tor /+ 

In this appendix we give the reducible contribution A of each diagram for the 
Kis form factor /+. The upper index of A refers to a specific diagram in fig. 2. 

The basic 1-loop functions occuring here (and in the other form factor /_) 
are Aq and Bq, and tensors and mass derivatives thereof, cf. appendix 
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II 

+ 

o' 

1 

(110) 

a( 1 “)/+ = 

{^ 21 ( 9 ^;™^.”!^) + 

(111) 

A(i'>)/+ = 

|3Ao(m^) + TAoCm^) + 5Ao(mJ)} 

(112) 

A<i‘=)/+ = 

+ |4L4(m^ + 2m^) + 2L^{-m?^ + m\) + g^Lgl 

(113) 

A(=“)/+ = 

■+ | 2 B 3 i(( 3 r^;m^,m^) - m^) 

(114) 


-2^31 mj) [8Li + 4 L 2 + L3]{m'^ — m^) 

+B2i(q^;Tn^, m^) [L3(m^ — + 3q^) — 6L4(mJ + 2m^) — 2L5(mJ + 5m^)] 

—-S 2 i(g'^; m^) [ 8 -Li(m^ — — q^) + 4L2Cm^ — + ij^) 


+L 3 (m^ — — 3g^) + 2L4(7m^ + lOm^) + 6 L 5 (m^ + m^)] | 


a(='>)/+ = 

— |3B2i(9^; iKk’ [ “ 4B4(m^ + 2 m^) - 2L^(m^ + m^) - ij^Lg] 

+B 2 i(g^; m^) [ — 12 Z/ 4 (mJ + 2 m^) + 2 L 5 (mJ — 7m^) - 39 ^ 1 . 9 ] | 

(115) 

A(^‘=)/+ = 

4lo(™^) [- 

3F<‘ 1 

-48Bim^ — 14B3m^ + 4L4(14m^ + mj) + 13B5(2m^ + mj) + 39 ^ Lg] 

(116) 


+Ao('m^) [—192Lim^ — 24L2r?T,^ — dOL^m^ + 4L4(44m^ + 7m^) + 2L5(31m^ + 7m^) + 7q'^L9] 
+Ao(r?T,^) [—144Lim^ — 24L2m^ — 54L3m^ + 4L4(10m^ + 29m^) + 2L5(2m^ + 23m^) + 5g^Lg] 
-2 A 2 (m^) [24L2 + 5 L 3 ] - 12 A 2 (m%) [4Li + I 8 L 2 + 5 L 3 ] - 6 [ 8 L 1 + 28 L 2 + 7 L 3 ]} 

/^(2d)f_^ ^ ^ 112 B 2 i(g^; m^) [L4(2m^ + mj) + Lsm^] (117) 

+4 B2i(q^; m%, m^) [3L4(2m^ + m^) + L5(4m^ — m^)] 

+ 12 B2i(q^; mj) [L 4 ( 2 m^ + mj) + LsmJ] 

+ 12 B2i(q^;Tn^, m%) [L 4 ( 2 m^ + mj) + L^m^] 

+ 12 B2i;i2(q^; ^k) [-^ 4 ( 2 m^ + mj) + L^rn]^ — 2LQ{2m^ + mj) — 2 L 8 m^] 

+ 12 -B 21 ; 12 ( 9 ^; [-^ 4 ( 2 m^ + m^) + Lsm^ — 2 L 6 ( 2 m^ + m^) — 2L3m'^] 

+ 12 B2i-,i2(q'^; rn^; m^) [L 4 ( 2 m^ + m^) + L^rn]^ — 2LQ{2m^^ + mj) — 2 L 8 r?T,^] 

+4 B 21 ; 12 (q^; [3L4m^(2m^ + m^) + L5m^(4m^ — m^) — 2Le{8m^ + 2m^m^ — m^) 

— 16L7(m^ — — 2L8(8m^ — 8 m^m^ + 3m^)]| 

I - 4Ao(m^) [3L4(2mK ++) + i5(4m^ - m^)] (118) 

— 28Ao{rri^j^) [L 4 ( 2 m^ + mj) + Lsm^] 

— 20j4o(m^) [L 4 ( 2 m^ + + L^rn^] 

-28Ao.2(m^) [L 4 ( 2 m^ + m^) + — 2LQ{2m\^ + m^) — 2 L 8 m^] 

— 20 Ao; 2 (?^^) [7;4(2m^ + m^) + L^rn^ — 2LQ{2m^ + — 2 L 8 m^] 

—4^0;2(?ti^) [3L4m^(2m^ + mj) + L5m^(4m^ — mj) — 2 L 6 ( 8 m^ + 2 m^mJ — m^) 

— 16L7(m^ — — 2L8(8m^ — 8 m^m^ + 3m^)]| 
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= + 2S2l(9^; rn^) B2i{q^m% , ml) 

(119) 





= {34lo(m^) B2i(q'^-,mp,m]i) + 3Ao(m^) B^iiq^, mp) 

(120) 


+9Ao(m%) B 2 i(q^-,m^,m%;) + 7Ao(m^) B 2 i{q^; m%,ml) 
+3Ao(m^) B2iiq^; trp, ni%:) + 5Ao(m^) B2i(<?^; 


a(3=)/+ 

- 24+4 {^B2i{q ■,m^,m^)[Ao{m^)+6Aoimj^)+Ao{m^)] 

(121) 


+B2i{q^; mp, m^) [33 Ao(m^) + 30 Ao(m^) + 9 Ao(m^)] 

-6 B2i-i2(q'^\ mp,-, mp) [Ao(m^) — 3 Ao(m^)] 

+3 B21;12(<?^; "Ik) + m^) 

+3 B21;12(<?^; mK) '4o(m^)(3m^ + m^) 



-2 B2i,i2{q^; m'\,; m^) [9m^ Ao(mJ) — 6(3m^ + mj) Ao(m^) + (16m^ - 

■ 7ml) Ao(m^)} 


1 1 

72 F4 1 

(122) 


57 Ao(m^) Ao{m'^) + 42 Ao(m^)^ + 41 Ao(m^) Ao(mJ) + 40 Ao(; 

2 n2 
^TT J 


-2 Ao; 2 (m^) Ao(m^) (16m^ - 7m^) + 12 Ao; 2 (m^) Ao(m^) (3m^ 
-18 Ao;2(m^) Ao(mJ) + 7 Ao.2{m\.) Ao(m^) (3m^ + mj) 

-10 Ao,2(m^) Ao(m^) + 30 Ao;2(m^) Ao(m^) m^j 


A (3-=)/^ 

- 18+4 S2l(<? +6Ao(m,,) B2 i(9 ■,m^,m^) 

(123) 


+24Ao(m^) B2 i(<?^; "1k) + 19Ao(m^) B 2 \(q^ \ m%, mj) 

+ 12Ao(mJ) B 2 i{q^; m^) + 20Ao(mJ) B 2 i(q^ \ mj)| 



Aoimlf+ 168 Ao{ml)Ao(7n^K) +'>-32 Ao{ml)Ao{ml) (124) 

+318 Ao(m%f + 386 Ao(m%) Ao(ml) + 175 Ao(m^)^} 


C Divergent parts of const ants 

In this appendix we list the divergent parts of all £*^®)-constants which occur 
in the meson vector form factors. They are derived from eqs. (fz^), (^) on 
the one hand and eqs. dzl), dH) on the other hand taking into consideration 
eq. dzl)- Similar relations for the electromagnetic form factors of K^, Kq, 
and the weak form factors f± for the rj —> decay are also taken into 

consideration [Ulil- 
Let 


Pi 

Pj = ^ 


(- 2 ) 


p: 


(- 1 ) 


+ pf + Oie) 


(125) 


e 

the Laurent expansion of the -parameters, cf. eq. (E^. For their divergent 
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parts 






£ 


( 126 ) 
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we find 


odiv 

P 22 


odiv I odiv 

P 22 P 23 


2i3^r-Pi^ = - 


2/3|r + /32t = 


/3it = 


div 


A 


odiv I odiv 

P 16 “t" P 18 


A: 


div 

17 


A: 


15 


1 


1 


(47r)2 [64(47r)2£2 


1 


1 


odiv 

P 26 


1 


_ - r(°) 
£ L96(47r)2 3 


" r I ^ r ^ r _i_ ^ r 

' 0^1 ^“^9 


i;-(0) 

4 ' 


(47r)2 2£"^3 

1 r 1 


+1^ 


7 


(47r)2 [ 32(47r)2£^ eLl92(47r) 

83 


1 


4 - 


-11 


1 


(47r)2 [ 96(47r)2£^ £Lll52(47r) 


V(o) V(o) 

— — ~JoQ 


- 2L^°^ + 


_ ^^(0) _ 1^(0) _ 


( 0 ) 


-( 0 ) 


If -7 1 r 35 

(47r)2|l28(47r)2e2 “^ £ L 1152(47r)2 3 


_i _r -I_ T 

+ 0^1 +6^ 


( 0 ) 
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( 0 ) 


-( 0 ) 


1 r -19 U 13 V(0) 

(47r)2 f 64(47r)2£2 £ L384(47r)2 3 ^ 3 ^ 


( 0 ) 


16 


_ - AO) _ “L 


( 0 ) 


1 


(47r)2 f 64(47r)2£2 

1 r 119 


(127) 

} 

(128) 

(129) 

(130) 

(131) 

(132) 

(133) 


+ 

£ L384(47r)2 3 ^ 


14. 


+ + 5L 


28 


( 0 ) 


-( 0 ) 


37 


- - QL)"' + &L 7 " - 2L 


18"® 
r -17 


( 0 ) 


( 0 ) 


( 0 ) 


(47r)2 f 128(47r)2£2 e L768(47r) 


173 


+ 7:L 


( 0 ) 


(134) 


+ 

3 ^ 


11 


^(0 )+Lr - 


( 0 ) 


13 

12 ' 


( 0 ) 


+ 12l1°^ + 3L 


( 0 ) 


7 


^ T 

(47r)2 f 192(47r)2£2 £ L324(47r)2 9 


( 0 ) 


A20 = 


+ 9^2 

1 f -11 


54 




( 0 ) 


1 AO) 
2^® 


1 


125 


(47r)2 \ 384(47r)2£2 e 1.2304(47r)2 9^^ 


(135) 


(136) 


7 


+ 


17 

los' 


(0) 


1 


AL® - 


17 


( 0 ) 


j^2“5 +4L7°^+2Z,g 


( 0 ) 
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where Lf''^ are the e°-coefficients of the -parameters, cf. eq. (il)- 

Note that the /3j do not contain any mass terms. Thus, the only source of 
masses in the Lagrangian dil) is the mass matrix x defined in 


D Electromagnetic pion form factor 

The 0{p^) calculation of the pion electromagnetic form factor is carried out in 
a manner completely analogous to that of the kaon semi-leptonic form factors. 
We quote the result of the contribution: 


^^+[£( 6 )] 


{/322 + 1323 } ( 137 ) 

2o2 

+ {-2m^/322 + 2 to ^/324 + m^/325 + {ml + 2m|-)/326 - 2 to ^/ 327 } ■ 


For the interpolation polynomials of the 0{p^) loop corrections (reducible plus 
irreducible diagrams) we obtain 

^ a:(0.0813-1-0.0619 a:-f 0.03873:^-1-0.0218 0.0068 cc^) 

(138) 


for X = lm\ G [—1, (2m^)^/m|.]) and in the decay region 

= - 0.0084-b 0.1076 a:-b 0.09560.0544 cc^-b 0.0233 (139) 

- 0.0047 X® -b i (0.0120 - 0.0998 x -b 0.2435 x^ -b 0.1848 x^ 

-b 0.1087x‘‘-b 0.0260 x^) 


for X = q^lm\ G [4m^/m|-, 1]). 

The separate reducible and irreducible 0{p^) loop contributions (modulo a 
quadratic polynomial) to the pion form factor are plotted in fig. 6. 

The arbitrary linear term can be htted by using the experimental pion charge 
radius 

= (0.439 ± 0.008) W . (140) 


The constant multiplying is obtained by using the curvature of the form factor 
at the origin. We chose a value of 3.4 GeV""', which represents an average of 
some typical extractions from experiment Q, Q, |21 . As fig. 7 demonstrates 


this choice of parameters yields a good description of the experimental form 
factor. 
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fig. 3. 



Figure 7\ The pion electromagnetic form factor in and 0(p^) chiral perturbation 

theory versus experiment. The data are from ref. p^], | |23| , and 
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